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ABSTRACT 



Recent observations have revealed that red, optically-passive spiral galaxies with 
little or no optical emission lines, harbour significant amounts of dust-obscured star 
formation. We propose that these observational results can be explained if the spa- 
tial distributions of the cold gas and star-forming regions in these spiral galaxies are 
significantly more compact than those in blue star-forming spirals. Our numerical sim- 
ulations show that if the sizes of star-forming regions in spiral galaxies with disk sizes 
of Rd are ~ 0.3i?d, such galaxies appear to have lower star formation rates as well as 
higher degrees of dust extinction. This is mainly because star formation in these spi- 
rals occurs only in the inner regions where both the gas densities and metallicities are 
higher, and hence the dust extinction is also significantly higher. We discuss whether 
star formation occurring preferentially in the inner regions of spirals is closely asso- 
ciated with the stripping of halo and disk gas via some sort of environmental effect. 
We suggest that the "outside-in truncation of star formation" is the key to a better 
understanding of apparently optically-passive spirals with dusty star-forming regions. 

Key words: stars:formation - galaxies:spiral - infrared: galaxies - galaxies:evolution 



1 INTRODUCTION 

Since the discovery of significant numbers of galaxies in dis- 
tant (z ~ 0.2-0.5) clusters with a spiral morphology but 
with no apparent on-going star formation based on the ab- 
sence of any emission lines in their optical spectra (e.g., 
Couch et al. 1994, 1998; Dressier et al. 1999; Poggianti et 
al. 1999), the origin of these so-called "optically-passive" 
spirals (or "k-type" spirals) has received considerable atten- 
tion both observationally and theoretically (e.g., Bekki et 
al. 2002; Goto et al. 2003; Yamauchi & Goto 2004; Moran 
et al. 2006; Masters et al. 2010). For example. Goto et al. 
(2003) found that such passive spirals are located anywhere 
between 1 — 10 virial radii from the centres of clusters and 
suggested that their formation is closely associated with 
cluster-related physical processes. In contrast, Masters et 
al. (2010) recently found that passive spirals exist preferen- 
tially in intermediate density regimes, and that there are no 
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obvious correlations between their physical properties and 
their environment. 

A further important and yet puzzling observational re- 
sult is that some of the passive spirals contain significant 
amounts of obscured star formation (e.g.. Wolf e al. 2005; 
Wilman et al. 2008; Wolf et al. 2009). The star formation 
rates in these cases are a factor of ~ 4 lower than those 
in blue spirals with the same mass. More specifically, the 
ratio of the star formation rate inferred from their infrared 
emission (SFRir) to that inferred from their UV emission 
(SFRuv) is typically a factor of ~ 3 larger than that for blue 
spirals (Wolf et al. 2009), implying that the passive spirals 
have a significantly higher ( x ~ 2) level of dust extinction. 
However, it remains unclear how and when such dust ex- 
tinction occurs in spiral galaxies when their star formation 
rates are significantly lower. 

The purpose of this paper is to show how optically- 
passive spirals with lower yet substantial star formation 
rates can have higher degrees of dust extinction, based on 
numerical simulations of star-forming disk galaxies. In par- 
ticular, we demonstrate that if the sizes of actively star- 
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"^ The initial gas mass fraction of a disk galaxy. 

The initial gas disk size of a disk galaxy measured in units of R^, where R^ (=17.5 kpc) is the initial (stellar) disk size of the galaxy. 
^ The initial mean surface gas density of a disk galaxy within ijg. 

The initial mean surface gas density of a disk galaxy within R^. 
^ The value of the extinction parameter a^ averaged over all time steps in a simulation. 
■< The value of the extinction parameter Ov averaged over the last 0.1 Gyr in a simulation. 
9 The value of the star formation rate averaged over all time steps in a simulation. 

The value of the star formation rate averaged over the last 0.1 Gyr in a simulation. 



forming regions [Rat) in galaxies with disk sizes of R^ are 
significantly more compact than i?d, then they will exhibit 
both lower star formation rates and heavier dust extinction. 
We assume that i?sf varies amongst spiral galaxies in the 
local and distant universe, and thus we treat R^i as a free 
parameter in this study. This assumption is consistent with 
the recent observations of Koopmann & Kenny (2004) which 
showed that isolated spirals as well as those in the Virgo 
cluster are extremely diverse in the radial distribution and 
extent of their Ha emission (from star- forming regions). We 
discuss how R^i can be more compact in some star-forming 
passive spirals in §4. 



2 THE MODEL 

We used the latest version of GRAPE (GRavity PipE, 
GRAPE- 7) - which is the special-purpose computer for grav- 
itational dynamics (Sugimoto et al. 1990) - in order to in- 
vestigate the chemodynamical evolution of star-forming disk 
galaxies. We have revised our original GRAPE-SPH code 
(Bekki 2009) for galaxy-scale hydrodynamical evolution so 
that we can investigate chemical evolution and star forma- 
tion processes of disk galaxies; the details of this new code 
will be described in future papers (e.g., Bekki et al. 2010). 
The masses of the disk, bulge ,and dark halo compo- 
nents of our model galaxy are represented by Md, Mb, and 
Afdm, respectively. The mass ratio of A/dm to Md was fixed 
at 16.7 for all of the present models so that the models 
can mimic the mass distribution of the Galaxy with the to- 
tal mass of ~ lO^^Mg (e.g., Evans & Wilkinson 2000). We 
adopted an NEW halo density distribution (Navarro, Erenk 
& White 1996) suggested from CDM simulations: 
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where r, po, and r^ are the spherical radius, the characteristic 
density of a dark halo, and the scale length of the halo, 
respectively. We adopted c = 7.8 (the ratio of rvir to rg, 
where rvir is the virial radius) that is reasonable for the 
total masses (~ IO'^^Mq) investigated in the present study. 
The radial (R) and vertical (Z) density profiles of 
the disk (with size Rd) were assumed to be proportional 
to exp{—R/Rdfi), with scale length Rd,o = 0.2_Rd, and 
sech^{Z/Zd,o), with scale length Zd,o = 0.04d in our units, 
respectively; both the stellar and gaseous disks follow this 
exponential distribution. In addition to the rotational veloc- 



ity caused both by the gravitational fields of the disk and 
dark halo components, the initial radial and azimuthal veloc- 
ity dispersions was assigned to the disk component according 
to the epicyclic theory with Toomre's parameter Q = 1.5. 

The mass ratio of the bulge to the disk (/b) and the 
scale length (.Rb.o) of the stellar bulge represented by the 
Hernquist profile were fixed at 0.167 and 0.04i?d, respec- 
tively, for all models, which is consistent with that of the 
bulge model of the Galaxy. In the present study, we describe 
only the results of the models with Md = 6 x 10^" M0 (thus 
Afb = 1-0 X 10^°Mq) and R^ = 17.5 kpc (thus _Rb,o = 0.7 
kpc). 

The cold interstellar medium (ISM) was distributed 
within Rg and modeled using SPH particles. The mass frac- 
tion of the isothermal ISM in the disk (/g) was assumed to 
be a free parameter. The initial temperature of the ISM was 
assumed to be 10* K. We adopted the same method as that 
used in Bekki & Chiba (2005) for determining the radial de- 
pendence of gas mass fraction in a disk for a given /g. We 
adopted the same chemical evolution model as those used in 
Bekki & Chiba (2005) and the chemical yield and the return 
parameter were set to 0.02 and 0.3, respectively. The star 
formation was assumed to follow the Schmidt law (Schmidt 
1959) with an exponent of 1.5. Kinetic energy of 10^^ erg 
per supernova is given to the ISM immediately after star 
formation occurs. 

Friel (1995) has derived the metallicity (Z) gradient of 
the Galactic stellar disk based on the ages and metallicities 
that are estimated for the Galactic open clusters. We there- 
fore allocated metallicity to each disk star according to its 
initial position as follows: 



[m/H]„ = [m/H],„ + adXi?, 



(2) 



where [m/H]^ q is the central metallicity. If we adopt plau- 
sible values of —0.091 for the slope a^ (Eriel 1995) and the 
central value of 0.48 for [m/HJ^ g, the mean metallicity of 
the disk is 0.0 in [Fe/H]. 

In order to more quantitatively estimate dust extinction 
around each individual star in star-forming disk galaxies, we 
introduced a dimensionless parameter, a^, which measures 
the degree of dust extinction for each new i-th stellar particle 
as follows: The dust extinction at wavelength A (Ax) around 
a star is described as follows (e.g., Spitzer 1978): 

AA = -2.51og-^ = 1.086A^dQccrd, (3) 

where _F^, fi/(0), A^d, Qc, and (Jd are the observed radia- 
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tive flux, the radiative flux in the absence of extinction, the 
column density (per cm^) along the line of sight, the di- 
mensionless extinction efficiency factor, and the geometrical 
cross section of a dust particle. 

We assumed here that aa and the dimensionless factor 
Qe are constant for all the models considered in this present 
study. Furthermore, previous models have shown that the 
dust mass of the Galaxy is linearly proportional to the total 
mass of the heavy elements (e.g., Dwek 1998). Therefore we 
took the degree of dust extinction (A^) around a new star to 
be proportional to A'd, which in turn is proportional to p^Z, 
where pg and Z are the 3D gas density and the gaseous 
metallicity around the star, respectively. Thus we defined 
the dust extinction parameter (a^.i) for each individual i-th 
new star as follows: 



Pg,iZi, 



(4) 



where pg,,i and Zi are the 3D gas density and metallicity 
around the star. We here consider that the present model 
enables us to discuss the importance of initial gaseous distri- 
butions in determining the mean dust extinction of a galaxy 
without using our previous fully consistent model consider- 
ing 3D dust distributions (Bekki & Shioya 2000). 

We mainly investigated the time evolution of the mean 
Ov of new stars and the star formation rate (bai) over a 
0.28 Gyr time interval in each simulation. We estimated the 
mean values of a^ and 6af in star-forming disk galaxies over 
the last 0.1 Gyr (av.n and 6sf,n, respectively) and adopted 
them as reasonable indicators of the amount of dust extinc- 
tion and global star formation for the galaxies. We consider 
that these mean values are better than the mean a^ and 
6af estimated for all time steps (ov and 6af, respectively), 
because strong starbursts can occur initially in the inner 
regions of the disks for most models, owing to very high 
gas densities there. Thus our estimates pertain to those 
times when star formation and chemical evolution proceeds 
steadily in the disks. 

The two key parameters in the present study are /g and 
Rg which control the sizes of the star- forming regions (Rsi)- 
Although we have run numerous models with different /g 
and _Rg values, we show mainly the results for two represen- 
tative models, the key parameter values for which are given 
in Table 1. This is mainly because these two comparative 
models show most clearly how the initial size of the gas disk 
is important in determining the mean star formation rate 
and the mean degree of dust extinction in a spiral galaxy. 
These two models are hereafter referred to as Model 1, which 
refers to a passive spiral with a lower star formation rate and 
higher level of dust extinction, and Model 2, which refers to 
a blue spiral with a higher star formation rate and lower level 
of dust extinction. We also briefly describe the dependences 
of mean a^ on model parameters (/g and 7?g). The mass and 
scale resolutions of the present simulations are 3 x IO^Mq 
and 193 pc, respectively, so that we can estimate Ov for local 
gaseous regions (~ 100 pc). 

Finally, we note that we are unable to discuss whether 
the simulated disks exhibit k-type spectra, because the 
present new chemodynamical simulation does not output 
spectrophotometric information. We will address this im- 
portant point in our future papers using chemodynamical 
simulations with spectroscopic synthesis code like those in 
our previous work (Bekki et al. 2001). 
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Figure 1. The distribution of new stars with ages less than 0.28 
Gyr projected onto the x-y plane (i.e., the disk plane) for the two 
representative models. Model 1 (left) and Model 2 (right). 




Figure 2. Time evolution of star formation rates (SFRs, upper) 
and mean a^ (extinction parameter, lower) for 0.28 Gyr in Model 
1 (filled black) and 2 (open green). 



3 RESULTS 

Figure 1 shows that owing to the initially different gaseous 
distributions, the distributions of newly-formed stars are sig- 
niflcantly different between the two models. Model 1 and 2. 
In Model 1, star formation can only proceed in the inner 
regions where both gas densities and metallicities are high. 
Owing to a larger gas mass fraction, new stars can be formed 
across the entire region of the disk in Model 2; star forma- 
tion can occur not only in the inner regions with higher 
gas densities and metallicities, but also in the outer regions, 
even though the gas densities and metallicities are lower. 
The mean star formation rate in Model 2 (lO.9M0yr~^ for 
the last 0.1 Gyr) becomes significantly higher than that in 
Model 1 (3.4M0yr~'^) owing to the initially larger /g. The 
spiral-like structure delineated by the very young stars with 
ages less than 0.28 Gyr in Model 2, suggests that well-defined 
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Figure 3. Plots of new stellar particles onto the Z — pg plane for 
Model 1 (left) and 2 (right). Here the metallicity (Z) and the local 
3D gas density (pg) around each individual new stellar particle is 
plotted. 
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Figure 4. The dependences of av,n (mean dust extinction for 
new stars with ages less than 0.1 Gyr) on Sg,g (initial mean gas 
density within Rg) for 8 different models. The models with Rg = 
O.SiJd and i?,g = R^ are shown in red and blue, respectively. 



large spiral structures can be more clearly seen in disks with 
globally active star-forming regions. 

In Figure 2 the time evolution of the mean star forma- 
tion rate and extinction for Models 1 and 2 are shown. Here 
it can be clearly seen that the total star formation rate fi- 
nally becomes as low as ~ 3M0yr~^ in Model 1, whereas 
its mean Ov becomes a factor of ~ 3 higher than that of 
Model 2 at T = 0.28 Gyr. The origin of this higher mean 
ttv in the disk of Model 1 can be explained via reference to 
Figure 3. This plots the 3D gas densities (pg) and metal- 
licities (Z) around the new stars formed within 0.28 Gyr in 
Models 1 and 2. Here it can be clearly seen that these two 
quantities are, on average, systematically higher in Model 
1 than in Model 2. This is due to the fact that in Model 

1 star formation occurs preferentially in the inner regions 
with higher gas densities and metallicities, whereas in Model 

2 it occurs even in the lower-density and lower- metallicity 
outer regions of the disk. These results suggest that the ob- 
served higher dust extinction in passive spirals with lower 
star formation rates might be closely associated with more 
centrally-concentrated gas distributions within them. 



In the present study, Eg,g depending both on /g and Rg 
is a key parameter that can determine av.n- Figure 4 shows 
that the mean a^ for the very young stars formed within the 
last 0.1 Gyr {av,n) depends on the initial mean gas densities 
within Rg (Sg,g) in such a way that av,n is higher in models 
with higher Sg^g. This figure also shows that models with 
Rg = Rd can show higher av,n, if /g ^ 0.3 (corresponding 
to Sg,g ^ 18.71M0pc^). Furthermore, Figure 4 shows that 
av,n can be quite low in models with Rg — 0.37?d if ./g is 
lower (^ 0.02): the star formation rates are also quite low in 
these models (< IMoyr"^). This result suggests that passive 
spiral galaxies need to have a certain minimum amount of 
gas centrally concentrated in their disks if they are to show 
both lower yet substantial star formation rates and higher 
dust extinction relative to normal blue spirals. 

We have shown that if star-forming regions are very 
strongly concentrated in the inner regions of spiral galaxy 
disks {Rg ^ 0.2Rd), then they appear to have rather low 
star formation rates (< IMQyr"^) and high levels of dust 
extinction (av,n > 0.15). For example, the model with /g — 
0.007 and Rg = 0.17?d shows bs!,n = O.85M0yr"^ and flv.n = 
0.30. This result implies that if most of the gas in gas-poor 
disk galaxies (with /g < 0.01) can be fueled to the nuclear 
regions and consumed rapidly there owing to some physical 
mechanism (e.g., galaxy-galaxy interaction), then such disk 
galaxies can be identified as passive spirals with nuclear star 
formation with higher degrees of dust extinction. 



4 DISCUSSION AND CONCLUSIONS 

If the scenario presented here (preferred star formation in 
inner regions of galaxies) is correct, this begs the question 
as to how the gas within disk galaxies might be truncated 
in this way in the course of their evolution. Previous nu- 
merical simulations show that ram pressure stripping by the 
hot intracluster medium can remove, quite efficiently, gas 
from the outer parts of spiral disks, so that their gas disk 
becomes much more compact than their stellar disk (e.g., 
Abadi et al. 1999; Kronberger et al. 2008). Furthermore, re- 
cent numerical simulations have shown that ram pressure 
stripping of halo gas, which is an important source of fuel 
for star formation in galactic disks, is more efficient in the 
outer parts of halos of disk galaxies in groups and clusters 
of galaxies (e.g., Bekki 2009). Thus it is possible that trun- 
cated gas disks can be formed as a result of halo and disk gas 
stripping, particularly in group and cluster environments. 

Recently, the Galaxy Zoo project has revealed a large 
optical bar fraction in red spirals at low redshift (70 ± 5% 
versus 27 ± 5% for blue spirals), and thus suggested that 
stellar bars are responsible for the truncation of star forma- 
tion in this subset of the spiral galaxy population (Masters 
et al. 2010). Tidal interactions between galaxies can trig- 
ger the formation of bars and consequently transfer rapidly 
disk gas into the inner regions of the galaxies (e.g., Noguchi 
1988). Therefore, the bars in disk galaxies can change the 
spatial distribution of gas such that the distribution can be 
much more centrally concentrated. Thus it could well also 
be possible that the origin of the proposed truncated gas 
disks has something to do with dynamical action of stellar 
bars in disk galaxies. 

An important and possibly testable prediction of the 
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scenario presented here is that any emission associated with 
the star formation in optically-passive red spirals should 
be much more compact that that associated with star for- 
mation in blue spiral galaxies. On the other hand, the star 
formation rate per unit area (i.e., star formation density 
measured in units of Mq yr~^ kpc^^) for these two differ- 
ent types of spirals should not be so different, because the 
star formation densities within the inner regions of passive 
spirals are expected to be as high as those in blue spirals. 
What is needed to test this is an emission line that traces 
star formation (and its rate) and is not heavily affected by 
dust extinction. Here, the Ha line probably holds the most 
promise in terms of being the optical line least affected by 
dust extinction and one that can be readily mapped spa- 
tially and out to high redshifts via high resolution imaging 
and integral field unit spectroscopy. 

As commented on above, a shortcoming of this study 
is its inability to show that disk galaxies with centrally- 
concentrated star formation exhibit passive k-type spectra. 
Although previous theoretical studies based on one-zone 
models and numerical simulations showed that e(a) and 
a-1-k/k+a spectra can be formed in dusty star-forming galax- 
ies (e.g., Shioya & Bekki 2000; Bekki et al. 2001; Shioya et 
al. 2002), they did not clearly show k-type spectra can be 
formed from dusty star-forming galaxies. Thus it is crucially 
important that as a next step we conduct further numeri- 
cal simulations that include spectrophotometric modeling 
which allow us to predict the spectroscopic signature associ- 
ated with star formation that proceeds in the inner regions 
of disk galaxies. 

If the star-forming regions of red passive spirals are as 
spatially extended as those in blue star-forming spirals, then 
it will be necessary to consider alternative scenarios to the 
one presented here. One such possibility worthy of brief men- 
tion here is that the upper-mass cutoff (m-upp) of the initial 
mass function (IMF) is significantly smaller (e.g., < 20Mq) 
in passive spirals. In this truncated IMF scenario, there are 
few or no massive O stars that can ionize the ISM (i.e., 
> 2OM0), such that (i) optical emission lines are very weak, 
and (ii) dust in the ISM can obscure star formation quite 
efficiently due to there being little destruction of dust by 
ionizing photons. This truncated IMF scenario has observa- 
tional support through being able to explain the UV and 
Hq properties of low surface brightness galaxies with low 
star formation rates (Meurer et al. 2009). More quantitative 
investigation based on numerical simulations of disk galaxy 
evolution with non-universal IMFs are required to test the 
viability of this scenario. 

Recent observational studies of distant galaxies based 
on Spitzer 24/im photometry and optical imaging by the 
Hubble Space Telescope have revealed the dusty nature of 
red galaxies and have provided new clues to the possible 
gradual truncation of galactic star formation in different en- 
vironments (e.g., Gallazzi et al. 2009; Woff et al. 2009). The 
present study suggests that truncation of star formation can 
occur more dramatically in the outer parts of disk galaxies, 
where environmental processes (e.g., tidal and ram pressure 
stripping) can be more effective. It also suggests that the 
possible inner dusty star-forming regions in passive spirals 
would be due to "outside-in truncation of star formation" in 
the course of disk galaxy evolution, in particular, in groups 
and clusters of galaxies. 



5 ACKNOWLEDGMENT 

We are grateful to the referee Christian Wolf for valuable 
comments, which contribute to improve the present paper. 
KB and WJC all acknowledge the financial support of the 
Australian Research Council throughout the course of this 
work. Numerical computations reported here were carried 
out both on the GRAPE system at the University of West- 
ern Australia and on those kindly made available by the 
Center for computational astrophysics (CfCA) of the Na- 
tional Astronomical Observatory of Japan. 



REFERENCES 

Abadi, M. C, Moore, B., Bower, R. G., 1999, MNRAS, 

308, 947 

Bekki, K., 2009, MNRAS, 399, 2221 

Bekki, K., Shioya, Y., 2000, ApJ, 542, 201 

Bekki, K., Chiba, M., 2005, MNRAS, 356, 680 

Bekki, K., Shioya, Y., Couch, W. J., 2001, ApJ, 547, L17 

Bekki, K., Shioya, Y., Couch, W. J., 2002, ApJ, 577, 651 

Bekki, K., Shioya, Y., Couch, W. J., Vazdekis, A., 2005, 

MNRAS, 359, 949 

Couch, W. J., Ellis, R. S., Sharpies, R. M., Small, I., 1994, 

ApJ, 430, 121 

Couch, W. J., Barger, A. J., Small, I., Ellis, R. S., Sharpies, 

R. M., 1998, ApJ, 497, 188 

Dressier, A., Small, I., Poggianti, B. M., Butcher, H., 

Couch, W. J., Ellis, R. S., Oemler, A., Jr., 1999, ApJS, 

122, 51 

Dwek, E., 1998, ApJ, 501, 643 

Evans N. W., Wilkinson M. I., 2000, MNRAS, 316, 929 

Friel, E. D. 1995, ARAA, 33, 381 

Gallazzi, A. et al. 2009, ApJ, 690, 1883 

Goto, T. et al., 2003, PASJ, 55, 757 

Koopmann, R. A., Kenney, J. D. P., 2004, ApJ, 613, 851 

Kronberger, T., Kapferer, W., Ferrari, C, Unterguggen- 

berger, S., Schindler, S., 2008, A&A, 481, 337 

Masters, K. et al., 2010, MNRAS, 404, 792 

Moran, S. M., Ellis, R. S., Treu, T., Salim, S., Rich, R. M., 

Smith, G. P., & Kneib, J- P., 2006, ApJ, 641, L97 

Navarro, J. F., Frenk, C. S., White, S. D. M., 1996, ApJ, 

462, 563 (NEW) 

Meurer, G. R. et al., 2009, ApJ, 695, 765 

Noguchi, M. 1988, A&A, 203, 259 

Poggianti, B. M., Small, I., Dressier, A., Couch, W. J., 

Barger, J., Butcher, H., Ellis, E. S., Oemler, A., Jr., 1999, 

ApJ, 518, 576 

Schmidt, M. 1959, ApJ, 129, 243 

Shioya, Y., Bekki, K., 2000, ApJ, 539, L29 

Shioya, Y., Bekki, K., Couch, W. J. 2001, ApJ, 2001, 558, 

42 

Spitzer, L., 1978, Physical processes in the interstellar 

medium. New York Wiley-Interscience. 

Sugimoto, D., Chikada, Y., Makino, J., Ito, T., Ebisuzaki, 

T., Umemura, M., 1990, Nat, 345, 33 

Wilman, D. J. et al. 2008, ApJ, 680, 1009 

Wolf, C, Gray, M. E., Meisenheimer, K., 2005, A&A, 443, 

435 

Wolf, C. et al., 2009, MNRAS, 393, 1302 

Yamauchi, C, Goto, T., 2004, 352, 815 



